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Charged lepton mixing as a consequence of neutrino mixing is studied for two generations e, ^ 
in the temperature regime <C T "C Mw in the early Universe. We state the general criteria 
for charged lepton mixing, critically reexamine aspects of neutrino equilibration and provide ar- 
guments to suggest that neutrinos may equilibrate as mass eigenstates in the temperature regime 
prior to flavor equalization. We assume this to be the case, and that neutrino mass eigenstates 
are in equilibrium with different chemical potentials. Charged lepton self-energies are obtained to 
leading order in the electromagnetic and weak interactions. The upper bounds on the neutrino 
asymmetry parameters from CMB and BBN without oscillations, combined with the fit to the solar 
\^ ' and KamLAND data for the neutrino mixing angle, suggest that for the two generation case there 

I is resonant charged lepton mixing in the temperature range T ~ 5 GeV. In this range the charged 

. lepton oscillation frequency is of the same order as the electromagnetic damping rate. 

(N ■ 
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Neutrinos play a fundamental role in cosmology and astrophysics [l*], and there is now indisputable experimental 
' confirmation that neutrinos are massive and that different flavors of neutrinos mix and oscillate 0. S S S S US j 
T-H , thus providing evidence for new physics beyond the Standard Model. Neutrino oscillations in extreme conditions of 
temperature and density are an important aspect of Big Bang Nucleosynthesis (BBN) , in the generation of the lepton 
! asymmetry in the early Universe [ll IslIM lol ITo| . and in the physics of core collapse supernovae[^ Hlj . 

An important aspect of neutrino oscillations is lepton number violation, leading to the suggestion that leptogenesis 
can be a main ingredient in an explanation of the cosmological baryon asymmetry |l2j |. Early studies of neutrino 
> propagation in hot and dense media focused on the neutrino dispersion relations and damping rates in the temperature 

■ regime relevant for stellar evolution or big bang riucleosynthesislH 0, 0|. This work has been extended to include 
, leptons, neutrinos and nucleons in the medium 'l'^. Matter effects of neutrino oscillations in the early universe were 
' investigated in I15I.I1(t|| and more recently a field theoretical description of mixing and oscillations in real time has been 

provided in ref . ■ While there is a large body of work on the study of neutrino mixing in hot and dense environments, 
much less attention has been given to the possibility of mixing and oscillation of charged leptons. Charged lepton 
Oh, number non-conserving processes, such as/i— >e7;/.t^3e mediated by massive mixed neutrinos have been studied 
ft) ■ in the vacuum in re fs.lia. Il^ I20II . For Dirac neutrinos the transition probabilities for these processes are suppressed 
^\ by a factor mi/Ml^^^^^ The WMAPpljl bound on the neutrino masses nia < leV yields typical branching 

■ ratios for these processes B < 10"''^ making them all but experimentally unobservable. In this article we explore the 
possibility of charged lepton mixing in the early Universe at high temperature and density. In section we discuss 
the general arguments for charged lepton mixing as a result of neutrino mixing and establish the necessary conditions 

■ for this mixing to be substantial. We suggest that large neutrino chemical potentials may lead to substantial charged 
M lepton mixing. 

Without oscillations BBN and CMB provide a stringent constraint on the neutrino chemical potentials [S^, 113 ^a, 



> 



with —0.01 < < 0.22 , |^^,t| < 2.6. Detailed studies |8l l24 125ll26ll27j show that oscillations and self-synchronization 
lead to flavor equilibration before BBN, beginning at a temperature T ~ 30 MeV'25'l with complete flavor equilibration 
among the chemical potentials at T ~ 2MeV 24, 25] . Thus prior to flavor equalization for T > SOMeV there could 
be large neutrino asymmetries consistent with the BBN and CMB bounds in the absence of oscillations. We study 
whether this possibility could lead to charged lepton mixing focusing on two flavors of Dirac neutrinos corresponding 
to e, /Lt and for simplicity in the temperature regime where both are ultrarelativistic, with <C T <C Mw- In section 
^we discuss general arguments within the realm of reliability of perturbation theory suggesting that the equilibrium 
state is described by a density matrix nearly diagonal in the mass basis. In this section we also discuss caveats and 
subtleties in the kinetic approach to neutrino equilibration in the literature and argue that results on the equilibrium 
state are in agreement with the interpretation of an equilibrium density matrix diagonal in the mass basis. Our 
main and only assumption is that for T > 30 MeV neutrinos are in equilibrium and the density matrix is nearly 
diagonal in the mass basis, with distribution functions of mass eigenstates that feature different and large chemical 
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potentials. While this is not the only, it is one possible scenario for substantial charged lepton mixing that can be 
explored systematically. In section Hill we explore charged lepton mixing in lowest order in perturbation theory as 
a consequence of large asymmetries in the equilibrium distribution functions of mass eigenstates. In this section we 
also critically discuss possible caveats and suggest a program to include non-perturbative corrections in a systematic 
expansion. In section ITVI we summarize the main aspects and results of the article. 



I. CHARGED LEPTON MIXING: THE GENERAL ARGUMENT 

Charged lepton mixing is a consequence of neutrino mixing in the charged current contribution to the charged lepton 
self energy. This can be seen as follows: consider the one-loop self energy for the charged leptons. The off-diagonal 
self-energy is depicted in fig. ^ for the case of electron-muon mixing. The internal line in fig. ^ (a) is a 
neutrino propagator off-diagonal in the flavor basis, which is non-vanishing if neutrinos mix. In Fermi's effective field 
theory obtained by integrating out the vector bosons the effective interaction that gives rise to charged lepton mixing 
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FIG. 1: Off diagonal charged lepton self energy: (a) one loop W-boson exchange, (b) self-energy in the effective Fermi theory. 

A simple Hartree-like factorization yields 

2Gp_ _ 

-y|-eL7p < VeLVf_,L > = EL (2) 

where the brackets stand for average in the density matrix of the system. Eqn. ^ gives the charged-lepton mixing 
part of the self energy as 

= 7'" < VeLVt,L > 7m ■ (3) 

The Fermi effective field theory contribution to the self-energy is depicted in fig. Q (b). 
The focus of this article is to study two aspects that emerge from this observation: 

• Mixing: The propagating modes in the medium are determined by the poles of the full propagator with a self- 
energy that includes radiative corrections in the medium. The full self-energy for the charged leptons is a 2 x 2 
matrix (in the simple case of two flavors), and eqn. ^ yields the off-diagonal matrix element in the flavor basis. 
This is precisely the main study in this article: we obtain the charged lepton propagator including radiative 
corrections in the medium up to one loop in the electromagnetic and weak interactions. Neutrino mixing leads to 
off diagonal components of the propagator in the charged lepton flavor basis. We flnd the dispersion relation of 
the true propagating modes in the medium by diagonalization of the full propagator including one loop radiative 
corrections. The true propagating modes in the medium are admixtures of electron and muon states: this is 
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precisely what we identify as mixing. The results given by equations H2I3|I state quite generally that electron 
and muon states are mixed whenever the neutrino propagator is off-diagonal in the flavor basis. We highlight 
that this is precisely the condition for flavor neutrino oscillations since the propagator yields the transition 
amplitude from an initial to a final state. Therefore we state quite generally that provided flavor neutrinos 
oscillate, namely if the neutrino propagator is off diagonal in the flavor basis, charged leptons associated with 
these flavor neutrinos will mix. The true propagating modes of charged leptons are linear superpositions of the 
charged leptons associated with the flavor neutrinos. We emphasize these statements because even though they 
are a straightforward consequence of flavor neutrino mixing, this precise point, and its consequences, have not 
been previously addressed in the literature. 

• Oscillations: Consider the decays W ^ i^eG or neutron beta decay n ^ peV in the medium. The electron 
produced in the medium at the decay vertex propagates as a linear combination of the true propagating modes 
in the medium, each with a different dispersion relation. Upon time evolution this linear superposition will 
have non-vanishing overlap with a muon state yielding a typical oscillation pattern. We study this oscillation 
between the electron and muon charged lepton by considering the evolution of an electron wave-packet produced 
locally at the decay vertex. These oscillations are akin to the typical oscillation between flavor neutrino states 
and are a consequence of the one loop radiative correction depicted in fig.JQl with neutrinos in the medium. The 
transition probability from an initial electron to a muon packet oscillates in time. While in the case of almost 
degenerate neutrinos oscillations are associated with macroscopic quantum coherence because the oscillation 
lengths are macroscopically large, this is not a necessary condition for oscillations, which occur whenever the 
initial state is a linear superposition of the propagating modes. The example of neutron beta decay gives a 
precise meaning to the statement of charged lepton mixing: in the decay of the neutron the charged lepton that 
is produced is identified with the electron. This is the initial state, which in a medium will propagate as a linear 
combination of the propagating modes with an oscillatory probability of finding a muon. These oscillations are 
fundamentally different from the space-time oscillations possibly associated with quantum entanglement and 
discussed in references j28j. 

Eqn. |j3J) generally states that there is charged lepton mixing when the density matrix is off diagonal in the fiavor 
basis. This is equivalent to the statement of neutrino mixing. A simple example of a density matrix off diagonal in 
the fiavor basis is p = |0„i >< 0,„| with |0,„ > being the vacuum state in absence of weak interactions but with a 
neutrino Hamiltonian with an off diagonal mass matrix in the fiavor basis. This is the interaction picture vacuum of 
the standard model augmented by a neutrino mass matrix with flavor mixing. In the two fiavor case with 



with vi^2{x, t) the fields associated with the mass eigenstates, 

< 0„j|i^e(f,i)l7^(f',t')|0m >= cos6lsin6'[< Qni\v2{x,t)v2{x' ,t' > - < 0,„|i/i(x',t)l7i(f' ,t')|Om >] ■ (5) 

If the propagators for the mass eigenstates only differ in the masses, this difference leads to a very small self-energy. 
In a medium the fiavor off diagonal expectation value Q could be enhanced by temperature and or density. Therefore 
the general criterion for substantial charged lepton mixing in a medium hinges on just one aspect: a large off diagonal 
matrix element < VeiT'tiL >■ One possible case for which this condition is fulfilled is if the density matrix is nearly 
diagonal in the mass basis with large chemical potentials for the different mass eigenstates. 

We emphasize that this is only one condition for substantial charged lepton mixing and by no means unique, the 
analysis above shows that the most general condition is simply that < VcLV^l > be large. 

In the general case a full solution of a kinetic equation should yield the value of < VeL^^iL >■ ^/ an equilibrium state 
of mixed neutrinos is described by a density matrix nearly diagonal in the mass basis with distribution functions for 
the different mass eigenstates with large and different chemical potentials, then simple expressions for the equilibrium 
propagators allow an assessment of the charged lepton mixing self energy. Can this be the case?. 



[x, t) = cos 9 ui [x, t) + sin 6 V2 (x, t) , j/^ (a;, t) = cos 9 V2 (x, t) 



sin 9 vi {x, t) 



(4) 



II. ON NEUTRINO EQUILIBRATION 



A. Equilibration in the mass basis 



A system is in equilibrium if 



(6) 
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where p is the density matrix of the system and H the total Hamiltonian H = Hq + Hint with Hq the Hamiltonian 
in the absence of weak interactions but with a mass matrix and Hint = Hj^c + Hcc- I^i the absence of weak 
interactions, an equihbrium density matrix po commutes with iJo, therefore it is diagonal in the mass basis. The 
equilibrium density matrix cannot have off-diagonal matrix elements in the mass basis because these oscillate in time. 
Of course without interactions the system will not reach an equilibrium state, however, as is the usual assumption 
in statistical mechanics, provided the interactions arc sufficiently weak but lead to an equilibrium state, an almost 
free gas of particles in equilibrium is a suitable description, and the canonical density matrix in such case is of the 
form Po = e"^"/"^ (the grand canonical could also include a chemical potential for conserved quantities). Examples of 
this are abundant, a ubiquitous one is the cosmic microwave background radiation: the Planck distribution function 
describes free photons in equilibrium, although photons reach equilibrium by undergoing collisions with charged 
particles in a plasma with cross sections much larger than those of neutrinos. 

Consider how the density matrix is modified from the "free field" form by "switching on" the weak interactions 
in perturbation theory. A perturbative expansion in the interaction picture of Hq begins by writing the interaction 
vertices in terms of neutrino fields in the mass basis. Neglecting sterile neutrinos, neutral current interaction vertices 
are diagonal and only the charge current interactions induce off-diagonal correlations in the mass basis. Let us write 
the full density matrix as p — po + Sp where 6p has a perturbative expansion in the weak coupling. The equilibrium 
condition leads to the following identity 

[Sp, Ho] = - [po, Hint] - [Sp, Hint] (7) 

Taking matrix elements in the mass eigenstates of Hq the solution of eqn. Q for the matrix elements of Sp in the mass 
basis can be found in a perturbative expansion. The matrix elements of Sp may feature non- vanishing off diagonal 
correlations in the mass basis as a result of charged current vertices which mix different mass eigenstates. However, 
a perturbative solution for the matrix elements of eqn. |(7J) in the mass basis would at most result in off diagonal 
correlations which are perturbatively small. Namely, the equilibrium density matrix is nearly diagonal in the mass 
basis. 

Expanding the field operators associated with the mass eigenstates in terms of Fock creation and annihilation 
operators of mass eigenstates, the spatial Fourier transform of the field operators is given by 

i^,(fc,0) =^a,(fc,A)Z^,(fc,A) + 6|(-fc,A)V,(-fc,A) ; i = 1,2 (8) 

A 

where the spinors U, V are orthonormalized positive and negative energy solutions solutions of the Dirac equation with 
mass rrii. If the density matrix is diagonal in the mass basis, then {al{k) aj{k)) cx Sij and the distribution functions for 
the mass eigenstates are (aJ(A;) aiik)) ; (&J(fc) bi{k)) for neutrinos and antineutrinos of mass i respectively. Switching 
on the neutral current interaction which is diagonal in the mass basis (provided there are no sterile neutrinos) will 
lead to the equilibration of neutrinos and the equilibrium distribution functions will be the usual Fermi-Dirac with a 
possible chemical potential. The charged current interactions yield vertices that are off-diagonal in the mass basis and 
induce cross correlations of the form (aj aj) with i ^ j. In free field theory this equal time correlation function, if non- 
vanishing, oscillates with a time dependence e'*-'^''""*-' *, however, in equilibrium there cannot be a time dependence 
of these off diagonal correlations as the following argument shows 

{a\{t) a,{t)) = Tr [pe^"' a\{Q) a,{Qy-'"*) = Tr (pal(O) a,(0)) (9) 

where we used eqn. Either the charged current interactions that generate these off-diagonal correlations exactly 
cancel the free field time dependence for all values of momentum k or, more likely, they lead to the decay of these off 
diagonal correlations to asymptotically perturbatively small expectation values as expected from the general arguments 
following eqn. 

This observation leads to the conclusion that if the perturbative expansion is reliable, the weak interactions lead 
to an equilibrium state described by a density matrix which is nearly diagonal in the mass basis but for possible 
perturbatively small off-diagonal elements. In perturbation theory the equilibrium distribution functions are diagonal 
in the mass basis and may feature a chemical potential for each mass eigenstate. Weak interaction vertices involve 
the flavor fields, but these are linear combinations of the fields that create and annihilate mass eigenstates, the true 
in-out states. Consider a far off-equilibrium initial state with a population of vector bosons (or neutrons) and no 
neutrinos, the decay of the vector bosons (or neutrons) results in the creation of a linear superposition of mass 
eigenstates, which propagate independently after production. CoUisional processes via the weak interaction lead to 
the decoherence of the mass eigenstates and ultimately to a state of equilibrium in which equal time expectation values 
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in the density matrix cannot depend on time. Neutral and charged current interactions yield different relaxational 
dynamics: in the mass basis the neutral current interaction is diagonal and relaxation processes via neutral currents 
lead to equilibration in the mass basis. Charged currents feature both diagonal and off-diagonal contributions in the 
mass basis, the diagonal ones yield relaxation dynamics similar to the neutral current interaction. The off diagonal 
contributions induce correlations between different mass eigenstates, but also lead to the relaxation of these off diagonal 
correlations. These two types of processes leading to relaxation dynamics for diagonal and off-diagonal correlations in 
the mass basis are akin to the different processes that lead to the relaxation times Ti (diagonal) and T2 (transverse) in 
spin systems in nuclear magnetic resonance |32j. These concepts are manifest in Stodolsky's effective Bloch equation 
description of neutrino oscillations in a medium with a damping coefficient in the "transverse" direction |3Cl|| whose 
inverse is the equivalent of the T2 relaxation time in spin systems. 

By the above arguments this asymptotic equilibrium density matrix must be nearly diagonal in the mass basis at 
least within the realm of reliability of perturbation theory. Equilibrium correlation functions of operators at different 
times must be functions of the time difference. Of particular relevance to the discussion below is the flavor off diagonal 
propagator (t'e(/c, t)V^{k, t')). Writing the flavor fields as linear combinations of the fields 1^1,2 this correlation function 
in the equilibrium density matrix diagonal in the mass basis is to zeroth order in the perturbation 

{ve{k,t)Vf,{k,t')) = -coaOsine {vi{k,t~t')Vi{k,0)) - {v2{k,t~t')V2{k,Q)) . (10) 

Mixed correlators of vi^2 cannot be functions of the time difference because of the different masses lead to a dependence 
ont + t'. As a simple but relevant example, consider the density matrix for the vacuum 

p=|0><0| (11) 

where |0 > is the exact ground state of H. This state can be constructed systematically in perturbation theory from 
the ground state |0m > of the Hamiltonian Hq in absence of weak interactions, namely the interaction picture ground 
state in the basis of mass eigenstates, 

|0 >= |0m > + \nm > 1 ■ (12) 

where In™ > are Fock eigenstates of Hq ("mass eigenstates") with energy Writing the full density matrix as 
P = |0m >< 0,„| + Sp one can find Sp systematically in perturbation theory. The off-diagonal flavor propagator 

Se^,ik,t-t') = < 0\iy,{k,t)V^{k,t')\0 >=< Orn\iyeik,t)V^{k,t')\Orn > +0{g) + ■■■ 

= cos 6* sin 61 < 0„|t/2(fc, t)T72(fc, i')|Om > ~ < 0™|j/i(fc, f)l7i(fc, <')|0m > +0{g) + ---. (13) 

This propagator is the lowest order intermediate state in the process We W fj,, and is also the internal fermion 
line, along with W- vector boson exchange in the off-diagonal self-energy contribution for charged leptons, see flg.lQ). 
This simple example also leads to conclude that if there is an equilibrium state for which the equal time "distribution 
function" < i'e(fc)l7^(/c) >^ then the unequal time correlation function 

Se^,{k,t-t')=<iyeik,tpf,{k,t')>^cos0sm0 < V2{k,t)v2{k,t') > - < vi{k,t)Vi{k,t') > =^0 (14) 

As it will become clear below, this is the correlation function that describes the charged lepton mixing. These 
arguments rely on the validity of the perturbative expansion and require revision in the case where perturbation 
theory in the mass basis must be reassessed. In section C|l we discuss this possibility and propose a method to 
re-arrange the perturbative expansion. 



B. On the kinetic approach 



Early kinetic approaches to the dynamics of oscillatin g n eutrinos in thermal environments were proposed by 
Dolgov 

13 HI 113, Stodolsky30i and Manoharj^. Dolgov |22j introduced density matrices in flavor space, whereas 
Stodolsky[23 and Manoharj^] used a single particle density matrix of flavor states leading to Bloch-like equations and 
a similar description was studied in ref. [13 . Stodolsky argued that decoherence between fiavor states emerged from a 
"transverse" relaxation akin to the relaxation time T2 in nuclear magnetic resonance '3 2]. A Boltzmann equation for 
mixing and decoherence was established by Raffelt, Sigl and Stodolsky in terms of a "matrix of densities" in the 
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nonrelativistic domain instead of the density matrix. In this approach the field operators for different flavors were 
truncated to only the annihilation operators and obtained a Boltzmann equation in a perturbative expansion. A fully 
relativistic treatment was presented in ref. [s^l introducing "matrices of densities" defined by the expectation value of 
bilinears of creation and annihilation operators of flavor states. A quantum kinetic description of oscillating neutrinos 
was presented in ref. |35j| with an approach similar to those of refs.|.30, ,33J in terms of single particle flavor states. 

All of the approaches to the kinetic description of oscillating neutrinos in a medium in one way or another use the 
notion of flavor Fock states, either in terms of single particle flavor neutrino states or by expanding field operators in 
terms of creation and annihilation Fock operators for flavor states. However, flavor states and theprecise definition 
of Fock operators associated with these states are very subtle and controversial |3a 1371 IsM l39l liol |41 | . 

Precisely, what are flavor states? While in the literature there is no precise definition of a "flavor state", 
a proper definition of such state should begin by expanding the flavor field operator in terms of Fock creation and 
annihilation operators. In such an expansion the spatial Fourier transform of the flavor neutrino field operator, for 
example the electron neutrino at t = 0, is given by 

i/e(fc, 0) = ^ ae{k, s)Ue{k, s) + f3l{-k, s)Ve{-k, s) (15) 

s 

where the spinors ; Ve are orthonormalized positive and negative frequency solutions of a Dirac operator with some 
mass and define a basis. A Fock state of an electron neutrino can be defined by 

|i/e(fc,s)) =at(fc,s)|0™) (16) 

where |0m) is the vacuum of the non-interacting theory, namely the vacuum of mass eigenstates. However, the 
expansion (|15(l requires a definite basis corresponding to a definite choice of the Dirac spinors U,V, which can be 
chosen to be solutions of a Dirac operator for any arbit rary mass. Each possible choice of mass gives a different 
definition of "particle". One possible choice is zero mass[3^, another choice would be the diagonal elements of the 
mass matrix in the flavor basis^^ o^' masses mi to the electron neutrino and m2 to the muon neutrino|33|- Any of 
these choices is just as good and physically motivated but obviously arbitrary. The creation and annihilation operators 
are extracted by proiection^42j, for example a\{k,s) = i'l(k,0)Ue{k, s). Writing the electron neutrino field operator 
as a linear combination of the field operators that create and annihilate mass eigenstates one finds 



al{k, s) = cos 9 



sin 9 



a{ (fc. A) Ul (fc, X)Ue (k, s) + bi {-k, A) V/ (-fc, X)Ue (k, s) 



a\{k, \)ul{k, \)Ue{k, s) + 62 (-fc, A) Vl{-k, \)Ue{k, s) 



(17) 



The transformation between the set of "flavor" operators and those that create and annihilate mass eigenstates is 
unitary and the scalar products of the spinors yield generalized Bogoliubov coefficients. It is clear that there is no 
single choice of spinor lAe that will make 

Ulik, \)Ue{k, s) = 1 ; Ul{k, X)Ue{k, s) = 1 

Vj(-fc, X)Ue{k, s) ; Vli-k, X)Ue{k, s) = (18) 

A surprising result of the above identification is that the annihilation operator ae{k, s) creates a linear combination 
of antineutrino mass eigenstates out of the vacuum |0m, > [37ll3^ I40ll4ll . 

This observation indicates that any choice of the solutions for the spinors Ue,Ve to define a flavor Fock creation 
operator leads to ^ cos0|i^i) + sui9\v2)- A possible definition of fiavor Fock states would be to define the 
flavor vacuum |0/ > as the state annihilated by the flavor annihilation operators (deflned for example for zero mass 
states) and to construct a Fock Hilbert space out of this vacuum by successive application of flavor Fock creation 
operators. However, while there is a formal unitary transformation that relates the flavor and mass Fock operators 
via the Bogoliubov coefficients, such transformation is not unitarily implementable in the infinite dimensional Hilbert 
space 37], in particular < 0/|0m 0. Finally one can simply define flavor states as 

\ve >= cos 6*11^1 > + sin 6*11/2 > • (19) 

However, because of the ambiguities with the deflnition of flavor Fock creation-annihilation operators discussed above 
these single particle states are indirectly related to the flavor fields Vi,,fj,{x,t) that enter in the standard model La- 
grangian. Furthermore, this single particle definition does not yield any information on a Fock representation of many 
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particle flavor states. A quantum statistical description of a neutrino gas is intrinsically a many body description, the 
total wave function of an n-fermion system must be completely antisymmetric under pairwise exchange. The second 
quantized Fock representation allows a systematic treatment of the many particle aspects, in particular in quantum 
statistical mechanics a distribution function is an expectation value of Fock number operator in the density matrix. 
Therefore simply defining flavor states as in ean. H19l) does not yield a complete information on the many particle 
nature of a neutrino gas. A systematic study of the many particle aspect of the dynamical evolution of a dense gas of 
flavor neutrinos with a physically motivated definition of flavor states even in the non-interacting theory was presented 
in ref. 42] wherein subtle but important effects associated with the non-trivial Bogoliubov coefficients in the dynamics 
were studied. Another subtlety emerges when a chemical potential is assigned to "flavor states" , a chemical potential 
is a thermodynamic variable conjugate to a conserved particle number. Even in the free theory, in absence of weak 
interactions, flavor number is not conserved if neutrinos mix and as a result a chemical potential for flavor neutrinos 
is not a well defined quantity even for free mixed neutrinos. Dynamical aspects associated with this issue were also 
studied in ref. [i^ . 

A counter argument to this critique would hinge on the fact that neutrino masses are small on the relevant energy 
scales and mass differences are even smaller, therefore one can approximately take all of the Dirac spinors to be 
practically massless. Of course if all spinors V for mass and flavor eigenstates are taken to be massless the overlaps 
(|18ll yield U^U ~ 1;U^V = and eqn. H19() becomes an identity. This point will become relevant below. 

While this approximation may be justified, it glosses over the main conceptual aspects and avoids the fundamental 
question of what precisely is a distribution function of fiavor states. Such function includes information over all scales, 
it yields the average occupation for all values of the momenta, not just the high energy limit. 

The main point of this discussion is that there are subtleties and caveats in the kinetic description based on "flavor 
states" or flavor matrix of densities, which involve Fock operators for flavor states. While these subtleties and caveats 
may not invalidate the broad aspects of the kinetic results, they cloud the interpretation of the equilibrium state of 
neutrinos. 

To highlight this point consider an equilibrium situation in which ne^i{k) =< v\{k)v^{k) >= 0, for a density matrix 
diagonal in the mass basis this means 



These conditions of "flavor equalization" are the same as those obtained in ref.[l5j for the equilibrium solution 
of the kinetic equations, although in that reference one active and one sterile neutrino were studied. The condition 
()20|l is consistent with identical chemical potentials for the mass eigenstates in the limit mi — m2 — 0. As discussed 
above, it is precisely taking mi = m2 = that yields the correspondence between the definition of the flavor states 
(|19|l and the relation between the flavor and mass eigenstates fields when all spinors are taken to be massless. This 
is also the approximation used in ref. |34| where flavor fields are expanded in the basis of massless spinors. These 
are precisely the approximations invoked in the kinetic approach and correspond to neglecting the neutrino masses. 
Restoring neutrino masses the off diagonal correlation in the mass basis would be Ue^ik) oc {rn\ — m\)lk'^. Thus 
an interpretation of the kinetic results is that the equilibrium state is described by a density matrix diagonal in the 
mass basis with equal chemical potential for the mass eigenstates with an off diagonal correlation rig^ cx (ra\ — m\) j k^ 
which is neglected in the kinetic approach. 

Thus the equilibriimr solution of the kinetic equation in ref. jlSjfsee eqn. (12) in this reference) can be interpreted as 
a confirmation of the statement of equilibration in the mass basis when the neutrino masses are neglected, although 
in the "flavor" formulation of the kinetic equations this information is not readily available. 

The discussion in the previous section based on general aspects of the full density matrix and a systematic per- 
turbative expansion avoids the caveats associated with the intrinsic ambiguities in the definition of flavor states and 
suggests that equilibration leads to a density matrix nearly diagonal in the mass basis. 

Quantum Zeno effect: References [s^, discuss the fascinating phenomenon of the Quantum Zeno effect 

or "Turing's paradox" 30J. In the case of neutrino mixing, this effect arises when the scattering rate is larger than 
the oscillation rate. Since neutrinos are produced in weak interaction vertices as "flavor eigenstates" when rapid 
collisions via the weak interactions which are diagonal in the flavor basis prevent oscillations, the states are effectively 
"frozen" in the flavor basis |30l|. This situation may be expected at high temperature. An order of magnitude estimate 
reveals that such a possibility is not available in the case under consideration, with a large difference in the neutrino 



cos6'sin0 < i>l{k)vi{k) > - < vl{k)v2{k) > = ^< v{{k)vi{k) >=< vl{k)v2{k) > 



(20) 



which in turn leads to the result 



Uee =< vl(k)Ve{k) > = < t'^(fc)i/p(fc) >= 71, 



(21) 
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asymmetries. The argument is the following: the oscillation frequency is given by 



where the matter potential 



V{k) KkGpT^L 



and L is the neutrino asymmetry difference between the two generations of neutrinos, 
possibility of large asymmetries, namely L ^ 1. The decay rates are of the order 



r ; 



(22) 



(23) 

Our study relies on the 



(24) 



The "quantum zeno effect" would operate provided F >> 
temperatures, the oscillation frequency 



f2. Even if V{k) >> Sm^ which can occur at high 



(25) 



and 1^ ~ GpT'^/L « 1 under the assumptions invoked in this article, namely: i) L ~ 1, ii) perturbation theory is 
valid in Fermi's effective field theory. A reversal of this bound would entail a breakdown of the perturbative expansion, 
and of the Fermi effective field theory. In perturbation theory the bound is even stronger if Sm^ >> V{k). Thus 
we conclude that, under the conditions studied in this article, the quantum Zeno effect is not effective in "freezing" 
the states in the flavor basis and that oscillations and relaxation may indeed result in a density matrix which is off 
diagonal in the flavor basis. 



C. Main assumptions 

After the above discussion on the general aspects of equilibration and the kinetic approach, we are in position to 
clearly state our main assumptions. These are the following: 

• i): for T » 1 MeV the electromagnetic and weak interaction rates ensure that leptons are in equilibrium in 
the early Universe, namely their distribution functions arc time independent. 

• ii): the results from the kinetic approach in refs.||25j indicate that for T >> 30 MeV neutrino oscillations are 
suppressed and flavor equilibration via oscillations is not operational. Therefore for T » 30 MeV neutrinos 
are in equilibrium, and there could be large asymmetries in the neutrino sector consistent with the BBN and 
CMB bounds in the absence of oscillations. 

• ill) : the arguments presented above lead us to assume that the equilibrium state of the neutrino gas is described 
by a density matrix which is nearly diagonal in the mass basis, and allow the distribution functions of mass 
eigenstates to feature different chemical potentials. As per the discussion above, this is consistent with the 
interpretation of the equilibrium state after equalization of the chemical potentials discussed after eqn. ipn|) . 

Equation H14|l entails that the equilibrium off-diagonal flavor propagator does not vanish. The numerical study 
in ref.[2§| shows that flavor equilibration with < 0.07 is established but for T ^ 2 MeV, well below the scale 
of interest in our study, clearly leaving open the possibility, which we assume here, of large asymmetries in the 
neutrino sector at a temperature much higher than that of flavor equalization. In summary: the combination 
of the general arguments suggesting that the equilibrium density matrix is nearly diagonal in the mass basis, 
at least within the framework of perturbation theory, along with the results of the kinetic approach lead us 
to assume that at high temperature T > 30 MeV neutrinos are in equilibrium, the density matrix is nearly 
diagonal in the mass basis and there could be large asymmetries for the different mass states consistent with 
the bounds in absence of oscillations. Equilibration of mass eigenstates implies that the neutrino propagators 
in the mass basis only depend on the time difference (translational invariance in time) and are determined by 
the equilibrium distribution functions, the off-diagonal flavor propagator is given by ean. H14|) . 

The dynamics that leads to equilibration and the mechanism by which substantial chemical potentials emerge is 
of course very important and require a much deeper and detailed investigation as well as an understanding of initial 
conditions for lepton asymmetries. A consistent study should address the subtleties and caveats associated with the 



9 



notion of flavor states or Fock operators. Such program is certainly beyond the scope and focus of this article. Here we 
study the consequences of this assumption in a perturbative expansion in the mass basis up to one loop order, namely 
up to leading order in the weak and electromagnetic interactions. In order to explore the main possible consequences 
of neutrino equilibration with large asymmetries within a simpler setting, we focus on the regime of temperature much 
larger than the charged Icpton masses. 

We emphasize, that our main observation, namely that charged leptons mix if neutrinos mix is one of principle and 
of a general nature. Our purpose is to study the potentially novel broad aspects of charged lepton mixing under these 
circumstances in the simplest scenario, postponing a more detailed analysis to a future study. 



III. EXPLORING THE CONSEQUENCES 



In perturbation theory in the electroweak interactions, the self-energy is computed in the basis of mass eigenstates, 
hence it is off-diagonal in the flavor basis if neutrinos are in equilibrium as mass eigenstates in the absence of 
oscillations. Including the electromagnetic self-energy, and introducing the charged lepton spinor fields with 
a , (3 = e, li the effective Dirac equation in the medium for the space-time Fourier transforms of these fields is the 
following [iTj 



^7'^ w — 7 • fc^ 



'q/3 



M^p + T.Zp{u, k) + J:^^{u;, k) + k)L + k) l\ i>p{oj, k) = 0, 

where M — diag(Me, M^t) is the charged lepton mass matrix and L = (1 — 7^)/2. 



(26) 



A. Electromagnetic Self-Energy 



The leading electromagnetic contribution to the charged lepton self-energy E^™(ijj, fc) for temperatures much larger 
than the lepton masses is dominated by one- photon exchange in the hard-thermal loop approximation [4^1 fi^]. As 
discussed below, the temperature region of interest for substantial charged-lepton mixing is T ~ Gev, therefore we 
will neglect corrections of order Mg /T^; M^/T^ ^ 10~^ (which already multiply one power of a) to leading order. 

Quark-lepton chemical equilibrium may lead to charged lepton chemical potentials as large as those for neutrinos, 
therefore we allow for arbitrary charged lepton chemical potentials with the possibility that fie,fi/T ^ 0(1). The 
self-energy E^™(a;,fc) is diagonal in flavor space I]'^™(ti;,fc) — diag (Sg™((jj, fc), I]^™(aj, fc)) . The matrix elements are 
given by if = e,ii)^^ 



x^Hfc 



+ 7- 



1 



with 



(27) 



(28) 



The function 



4 f 



xdx 



s^-«/ + 1 



1 



4/ - ^ 



(29) 



is monotonically increasing with —0.196 < D{^) < 0.4, for ~ I < £, < 1- 

The leading order damping rate F in Y."" emerges from threshold infrared divergences associated with the emission 
and absorption of soft, transverse magnetostatic photons 47, 42] and is insensitive to the masses for temperature 
larger than the mass of the charged leptons. In an abelian plas ma transverse photons do not acquire a magnetic 
mass and are only screened dynamically by Landau damping j46l l47l |48| . Detailed work in high temperature QED 
plasmas j46l l47l liq reveals that the exchange of soft magnetostatic photons yields an anomalous damping of fermionic 
excitations, which can be accurately described by the damping rate given by 



aT In f ^) 



(30) 



The plasma frequency LUp is determined by the photon polarization fermion loop with all relativistic fermionic species: 
for T ~ GeV with electron, muon and three light quark degrees of freedom we find to leading logarithimic order 



F = aTln ( — ) ~ 3.3 aT. 



(31) 
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w 




a a 

FIG. 2: Charged current self energy Ea(w, k) with a = 1, 2 corresponding to mass eigenstate neutrinos in the fermion line. The 
external lines correspond to e, charged leptons with frequency u and momentum k. 



A collisional contribution to the charged lepton damping rate F is of order a'^TWc\ and will be neglected to leading 
order in a. 



B. Charged and neutral currents self-energy 

As depicted in fig. Q to lowest order in perturbation theory in the interaction picture of Hq, there is a flavor off- 
diagonal contribution to the charged lepton self energy, Se,/i given by a W-boson exchange and an internal off-diagonal 
fermion propagator < i7e{x,t)i/fx{x' ,t') >. In the mass basis this propagator is given by eqn. 1)14(1 and the general 
form of the corresponding self energy contribution is depicted in fig|21 where the internal fermion line corresponds to 
a mass eigenstate neutrino Va in equilibrium with chemical potential iJ,a = T ^a- 

We focus on the temperature regime T <C Mw and obtain the charged current contribution to the self-energy 
up to leading order in a local expansion in the frequency and momentum of the external leptons neglecting terms 
proportional to nia/Mw < 10~^^. The general expressions for charged and neutral current self-energies are given in 
references (ij, [iTj where we refer the reader for more details. 

Denoting by T,^'-^ {lj, k) the one-loop charged current self-energy with internal neutrino line corresponding to a mass 
eigenstate Va displayed in figEl the matrix S^_'^(a;, k) in eqn. has the following entries 

^^^^ = S^Si + C^Sa 

= S;:,? = -C5(Si-E2) (32) 

consistently with a perturbative calculation in the mass basis to lowest order in the weak interactions C = cos 9 ; S — 
sinO and 9 is the neutrino vacuum mixing angle. A fit to the solar and KamLAND^^ data yields tan^ 9 ~ 0.40. 

The off diagonal element in eqn. (|32|) is responsible for charged lepton mixing. The flavor off-diagonal 

propagator in this self-energy is precisely given by eqn. H14|l . The form of this off-diagonal self-energy makes clear 
that there is mixing provided 9 ^ 0,n/2 and Si — S2 7^ 0. The calculation of the self-energy S^J^ in the vacuum 
is standard: it is performed in the interaction picture of the true basis of in-out states, these are mass eigenstates. 
In this case the difference of the self-energies is determined solely by the neutrino mass difference, therefore in the 
vacuum oc G^Am^ 10~^^ and the mixing between charged leptons is negligible. 

The main point of our study is that in the medium in equilibrium with large neutrino asymmetries for the mass 
eigenstates, charged lepton mixing may be substantial. The propagating modes in the medium are determined by the 
poles of the exact propagator. An off diagonal self-energy entails that the charged lepton propagating modes are 
admixtures of electron and muon degrees of freedom. We now study this possibility in detail when the temperature 
is much larger than the lepton masses. We focus on this case for simplicity in order to extract the main features of 
the phenomenon and to highlight the main steps in the calculation. 

We are only interested in the real part of the self-ener gies since at this order the imaginary part vanishes on the 
charged lepton mass shells. From the results obtained in[T3, for any loop with a lepton with mass m in the high 
temperature limit T 3> m we obtain 



i?eE(w,fc) =7° a"{uj,k)-^-ka\uj,k). 



(33) 
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For Mw.z > T, t^, fc we find 



1 TTT^gCpkr^. 



where g is the appropriate factor for charged or neutral currents, L is the asymmetry for the corresponding lepton 
and Mb = Mw,z for charged or neutral currents, = 2 C(3) T^/tt^, and 



120 



77r4 L e^~« + 1 



dx . (36) 



The neutrino asymmetries are given by La — (7r^/12C(3))(^a + Ca/7r^). In the absence of oscillations the combined 
analysis from CMB and BBN yield an upper bound on the asymmetry parameters |^e| < 0.1, |^^| ^ 1 for flavor 
neutrinos \l'A which we assume to imply a similar bound on the asymmetries for the mass eigenstates, ^a- The 
validity of this assumption in free field theory is confirmed by the analysis in ref . [4^ • From H32|l and the above results 
the following general form for the charged and neutral current self-energies is obtained 



i?eE(^,fc) = i|^[7°AH-7-kB(fc)] (37) 
where A(ijj, fc) and ]B(ijj, fc) are 2x2 matrices in the charged lepton flavor basis given by 

- ( tt! tt) ) • »<^-' ^ ( b:S £S ) ■ 

where the matrix elements are given by 



and 



L+ + cos26'L_ 



i+ - cos26'i_ - 



(/+ +cos26l/_) 



90 C(3) Ml, 



^f„^(^) = sin 20 



90C(3)M2, 



(/+ +cos26'/_) 



270 C(3) Ml 

— ^(/+ -cos26'/_ 

270 C(3) M2,^ + 

77^^^ kT . 

— TT- sm 2t^ i_ , 

270 C(3) 



(39) 
(40) 
(41) 

(42) 
(43) 
(44) 



where we have introduced L± — Li ± L2, I± = h ± I2. The neutral current contributions are flavor diagonal and 
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given by 



-(1-4 sin^l 
-(1-4 sin^ ( 
-4 sin'' 
-4 sin'' 6',, 



" 90 C(3) 
^ 90 C(3) 



90 C(3) 
90 C(3) 



and 



77r4 /fc T 



270 C(3) 
77r4 A: T 



270 C(3) 



(1 — 4 sm Oyj) 



cos^ 




^ f 


cos^ 




/ 




8 

+ 3 


sin^ei^ ^g)Lf 
f 




8 

+ 3 


sin^ 6^ ^ g]; L/ 
/ 


COS 






cos^ 






4 








(1-4 sin^ e^) ' 



(45) 
(46) 
(47) 
(48) 

(49) 
(50) 
(51) 



where , Lf are the vector coupUng and asymmetry of fermion species /, 9^ is the Weinberg angle and £e^^ are 
the charged lepton asymmetries. The non-vanishing ofT-diagonal matrix elements A'^^ , lead to charged lepton 

mixing and oscillations. It is convenient to combine the charged leptons into a doublet of Dirac fields 



(52) 



In the chiral representation the left and right handed components of the Dirac doublet are written as linear combina- 
tions of Weyl spinors w''*^ eigenstates of the helicity operator a ■ k with eigenvalues h = ±1, as follows [T^ 



h=±l 







The left handed doublet 



(^(")(u;,fc) 



E 

h=±l 







(53) 



(54) 



obeys the following effective Dirac equation in the medium to leading order in a, GfUj 

3 Gp 



[(u; + ^^)2-fc2] 1 



2V2 



2wA^^-2fc]B^^ + (t^-M)(A + /iB) -M^ \ (p^'''> {lu , k) ^ , (55) 



where = diag (Afg -I- 2ml, ^fl + where m^^^ are given by eqn. ipSjl and to avoid cluttering of notation 

A , B are the sums of the charged and neutral current contributions. To leading order the right handed doublet is 
determined by the relation^ 



(56) 



The propagating modes in the medium are found by diagonalization of the above Dirac equation. Let us introduce a 
doublet of collective modes in the medium 



(57) 
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related to the flavor doublet (p^^^ [uj, k) by a unitary transformation Um^ with 



jj{h) ^ 



COS 6; 



sin 



- sin 9m 



m "^uo Ufa 




(58) 



(59) 



and a similar transformation for the right handed doublet ^'-''^(o;, fc), where the mixing angle 0m •* depend on h, k and 
Lo. The eigenvalue equation in diagonal form is given by 



I (a; + iTf -k^ + l Shico, k) - ^ [Ml + + 2 + 2 ) + 1 n^,{u;, k) J 
where Sh{i-o, k), Ah{io, k) and ^lh{Lo, k) are respectively given by 



X^''\oj,k) = Q 



(60) 



Sh{uj,k) = 



(Cc> - hk) [A^^iio) + Aee(c^) + h{ B,,{k) + B^^(fc))] 



2a.(<,^ + Af,^)-2fc(Sil'^ + i?,^^) , 



(61) 



3G 



2a.«^ - Al^) - 2fc(<^ - Bg^) , 



r!;,(cc;,fc) = 5M2-Aft(cj,fc) +[2{Lo^hk){A,^, + hBe^)] 



(62) 
(63) 



where iJM^ = M,^ 



+ 2m^ — 2mg. The mixing angle in the medium is determined by the relations 
sin20W =-^"''^^^^^^ + ^-^^^^ ~ 



Vlh{uj,k) 



(64) 



where w must be replaced by the solution of the eigenvalue equation (|60|l for each collective mode. A remarkable 
convergence of scales emerges for T ^ 5GeV: if the neutrino asymmetry \L_\ ~ 1 then for nearly thermalized 
relativistic charged leptons with uj ^ —hk with k ^ T, all of the terms in the expression for ri?i(w, k) are of the same 
order, namely, \Ah{k, k)\ ~ |rAep| ~ SM^. For relativistic leptons |w| can be replaced by k in the arguments of the 
functions A, B to leading order in a,GF- 

For <C w, fc,T <C Mw we find that the leading contribution to is given by 



Ah(uj,k) ~ 1.210" 



GeV 



hk 



2T 



L_ COS26I-0.29 



6.3w — hk 
7!3T 




(GeV^) , 



(65) 



where we have used the value sin 9^ = 0.23. A resonance in the mixing angle occurs for Ah = SM'^. The typical 
momentum of a lepton in the plasma is fc ^ T, therefore in the temperature regime T <C Mw wherein our calculation is 
reliable, a resonance is available uje,fj, (k) ^ —hk —hT if the neutrino asymmetry is close to the upper bound. Ta king 
the values for inferred from the upper bounds from combined CMB and BBN data in absence of oscillations [23. |2^ 
and the fit tan^ 9 ~ 0.40 from the combined solar and KamLAND data 49] suggest the upper bound |L_ 008 26*1 ~ 1. 
With the asymmetry parameters for the charged leptons |^/| smaller than or of the same order of |^a|) resonant mixing 
may occur in the temperature range T ~ 5 GeV. Even when the asymmetries from charged leptons do not allow for 
a resonance or at lower temperature, it is clear that at high temperature and for large neutrino asymmetries such 
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that L_ ~ 1 there is a large mixing angle because of the convergence of scales. Hence at high temperature and large 
differences in the chemical potential for mass eigenstates, the propagating charged lepton collective excitations in 
the medium will be large admixtures of e ; /i states. Consider a slightly off-equilibrium disturbance in the medium 
corresponding to an initial state describing an inhomogeneous wave packet of electrons. The real time evolution of 
this state in the medium has to be studied as an initial value problem. Following the real time analysis presented in 
ref.^3|j fi^d that if an initial state describes a wave-packet of left handed electrons of helicity h, with amplitude 
?g (0; k); k 3> but no muons, the persistence and transition probabilities are given by 



Pe-.e{t;k) = rMk)\'e 



2 „-2rt 



Pe^^{t;k) = \l^A0;k)fe-^"sm\2et{k))si 



l-sin2(2<(fc))sin2 

n{k) 



n{k) 

4fc 



4fc 



; n{k)^nh{k,k). 



(66) 
(67) 



The exponential prefactor reveals the equilibration of the charged lepton distribution with the equilibration rate 
2r^^|45. It is also remarkable that T ^ Vt{k)/k in the temperature and energy regime of relevance for the resonance 
k T ^ b GeV. Therefore we conclude that during the equilibration time scale of charged leptons, there is a 
substantial transition probability. CoUisional contributions are of order a^T or Gj^T^ for electromagnetic or weak 
interaction processes leading to collisional relaxation time scales far larger than the oscillation scale for T 5 GeV. 

In the radiation dominated phase for M\y ^ T as discussed here, we find that the ratio of the oscillation to the 
expansion time scale Htqsc ^ 10~^^(T/GeV)'^ <C 1, namely oscillation of mixed charged leptons occur on time scales 
much shorter than the expansion scale and the cosmological expansion can be considered adiabatic. 



C. Remarks: beyond perturbation theory 

We have focused on the high temperature limit to provide a detailed calculation within a simpler scenario, to extract 
the main aspects of the phenomenon and to highlight the main steps of the calculation. However, it is clear that a 
similar calculation can be performed at much lower temperature and the point of principle is still valid under the 
assumption of an equilibrium density matrix diagonal in the mass basis: there could be substantial charged lepton 
mixing if there are large chemical potential differences between the distribution functions of mass eigenstates. As 
per the discussion above, this is not the only scenario that yields substantial charged lepton mixing, the general 
condition is that the off diagonal self energy Sg^ in eqn. be non-zero (and large). The off diagonal expectation 
value < VeLV^j^L > must in general be found from the equilibrium solution of a kinetic equation, but with a consistent 
treatment that avoids the caveats and subtleties discussed in section Ijll Bll . 

The arguments in favor of an equilibrium density matrix diagonal (or nearly so) in the mass basis, and the specific 
calculation described above relied on a perturbative expansion in the interaction picture of the unperturbed Hamilto- 
nian Hq which includes the neutrino mass matrix. There are possible caveats in the validity of perturbation theory, 
particularly in the case where medium effects lead to large corrections to the single particle states. For example, a 
large "index of refraction" arising from forward scattering with particles in the medium may lead to non-perturbative 
changes in the properties of the single particle basis. The lowest order contribution to the self-energy from forward 
scattering has been obtained in ref.[l2|, these are in general dependent on the energy of the neutrinos. Including these 
corrections in a perturbative approach entails summing the geometric Dyson series for the one-particle irreducible self 
energy in the neutrino propagators. This case is akin to the generation of a thermal mass from forward scattering in 
a scalar (f>'^ field theory at finite temperature ^465 , when this thermal mass is larger than the zero temperature mass 
there is a large modification in the propagating single particle modes in the medium. In the scalar field theory case a 
self-consistent re-arrangement of the perturbative expansion consists in adding the thermal mass to the unperturbed 
Hamiltonian and at the same time a counterterm in the interaction part. The free single particle propagators now 
include the thermal mass term, and in order to avoid double counting, the counterterm in the interaction Hamiltonian 
cancels the contributions that yield the thermal mass corrections systematically order by order in the perturbative 
expansion. We propose a similar str ategy to include the medium modifications to the propagating single particle 
modes in the medium. In references [13, Il7| it is found that the forward scattering contributions to the effective 
Hamiltonian in the medium are of the form 

= 7°A(fc)-7-fc]B(fc) (68) 
with A(fc); ]B(fc) momentum dependent matrices in flavor space, their explicit expressions are given in refs.|l.'ill7j. A 
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re-arrangement of the perturbative expansion results by writing 

H = Ho + H,„t (69) 

where Hq — Hq + 6H and Hint — Hint + He where the "counterterm" Hamiltonian He — ~6H systematically 
cancels the forward scattering corrections to the self-energies consistently in the perturbative expansion. The new 
"free" Hamiltonian Hq includes self-consistently the modifications to the propagating single particle states from the 
in-medium index of refraction. The field operators are now written in the basis of the solutions of the Dirac equation 
from the new Hamiltonian and finally the interaction is written in terms of these fields. Thus the perturbative 
expansion is re-organized in terms of the single particle propagating modes in the medium. The main complication 
in this program is that the mixing angles in the medium which determine the single particle propagating modes, 
are energy dependent, this introduces a non- locality in the interaction vertices which now become dependent on the 
energy and momentum flowing into the vertex. 

If there is substantial charged lepton mixing, such phenomenon in turn affects the index of refraction for neutrinos 
in the medium and possibly the equilibration dynamics of neutrinos. A self-consistent treatment of the charged lepton 
mixing and neutrino mixing and relaxation would be required to understand the dynamical aspects of neutrino and 
charge lepton equilibration. This task is beyond the goals and focus of this article. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this article we focused on studying the possibility of charged lepton mixing as a consequence of neutrino mixing 
at high temperature and density in the early Universe. There are three main points in this article: 

• (I) We establish that a general criterion for charged lepton mixing as a consequence of neutrino mixing is that 
there must be off diagonal correlation of flavor fields in the density matrix. We identified one possible case in 
which there could be charged lepton mixing: that the equilibrium density matrix be nearly diagonal in the mass 
basis. This is the case in the vacuum, but in this case the smallness of the neutrino masses entails that charged 
lepton mixing is negligible. We argued that this effect can be enhanced in a medium if the density matrix is 
nearly diagonal in the mass basis with large and different chemical potential for mass eigenstates. While this 
is not the only case in which charged lepton mixing can occur, it is one in which we can provide a definite 
calculation to assess charged lepton mixing. 

• (II) We have given general arguments to suggest that within the realm of validity of perturbation theory, 
the equilibrium density matrix must be nearly diagonal in the mass basis. We have critically re-examined 
the kinetic approach to neutrino mixing and relaxation in a medium at high temperature and density and 
highlighted several caveats and subtleties with flavor states, and or Fock operators associated with these states 
that cloud the interpretation of the density matrix. We argue that an equilibrium solution of the kinetic 
equations describing "flavor equalization" [l5| can be interpreted as a confirmation that the density matrix is 
nearly diagonal in the mass basis. This interpretation leads us to the main and only assumption, namely that 
before "flavor equalization" for T > 30 MeV neutrinos are in equilibrium, the density matrix is nearly diagonal in 
the mass basis but with distribution functions for mass eigenstates with large and different chemical potentials 
in agreement with the bounds from BBN and CMB in absence of oscillations. 

• (III) Under this assumption and the validity of perturbation theory we have provided a definite calculation of 
charged lepton mixing. While the general criterion for charged lepton mixing does not imply that this is the only 
case in which charged leptons mix, it is a scenario that allows a definite calculation to assess the phenomenon 
in a quantitative manner. 

In conclusion, under the assumption that the mass eigenstates of mixed neutrinos in the early Universe are in 
thermal equilibrium with different chemical potentials for T > > 30 MeV, before oscillations establish the equalization 
of flavor asymmetries |25j| neutrino mixing leads to charged leptons mixing with large mixing angles in the plasma. 
We explored this possibility by obtaining the leading order contributions to the charged lepton self energies in the 
high temperature limit. If the upper bounds on the neutrino asymmetry parameters from BBN and CMB without 
oscillations is assumed along with the fit for the vacuum mixing angle for two generations from the KamLAND 
data, we find that charged leptons mix resonantly in the temperature range T ~ 5 GeV in the early Universe. The 
electromagnetic damping rate is of the same order as the oscillation frequency in the energy and temperature regime 
relevant for the resonance suggesting a substantial transition probability during equilibration. The cosmological 
expansion scale is much larger than the time scale of charged lepton oscillations. Although we assumed the validity 



16 



of perturbation theory we recognized possible caveats in the high temperature hmit arising from potentially large 
corrections to the single-particle propagating modes from the in-medium index of refraction. We proposed a re- 
organization of the perturbative expansion that includes the correct single- particle propagators self-consistently. We 
have focused on the high temperature limit as a simpler scenario to assess charged lepton mixing, however, the 
calculation can be performed at lower temperatures with the corresponding technical complications associated with 
the lepton masses. While at much lower temperatures there is no resonant mixing of charged leptons, the results of 
the calculation establish a point of principle, namely that for large chemical potential differences in the distribution 
function of mass eigenstates of neutrinos, the charged lepton propagating modes in the medium will be admixtures 
of the electron and muon degrees of freedom with non- vanishing mixing angle. 

We believe that the phenomenon of charged lepton mixing in a medium warrants a deeper and thorough investi- 
gation. Our study also raises relevant questions on the kinetic approach: a consistent description of the kinetics of 
neutrino oscillation and relaxation avoiding the caveats associated with flavor states and or Fock operators associated 
with these states. Furthermore, substantial charged lepton mixing also suggests that a dynamical description should 
include self-consistently both neutrino and charged lepton mixing in a full non-equilibrium treatment. These aspects 
as well as possible consequences of charged lepton mixing for leptogenesis will be explored elsewhere. 
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